Comment on "Proposed central limit behavior in deterministic dynamical systems" 
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In a recent Brief Report [Phys. Rev. E 79 (2009) 057201], Grassberger re-investigates probability 
densities of sums of iterates of the logistic map near the critical point and claims that his simulation 
results are inconsistent with previous results obtained by us [U. Tirnakli et at, Phys. Rev. E 75 
(2007) 040106(R) and Phys Rev. E 79 (2009) 056209]. In this comment we point out several errors in 
Grassberger's paper. We clarify that Grassberger's numerical simulations were mainly performed in 
a parameter region that was explicitly excluded in our 2009 paper and that his number of iterations 
is insufflcient for the region chosen. We also show that, contrary to what is claimed by the author, 
(i) Levy distributions are irrelevant for this problem, and that (ii) the probability distributions of 
sums that focus on transients are unlikely to be universal. 

PACS numbers; 05.20.-y, 05.45.Ac, 05.45.Pq 

In a recent Brief Report [l|, Grassberger claims that his numerical re-investigation of the probability density of 
sums of iterates of the logistic map near to the critical point of period doubling accumulation is inconsistent with our 
results previously obtained in 0, Q ■ In 3we provided for the first time numerical evidence for the possible relevance 
of g-Gaussians for this problem, and in [3| a more detailed investigation was performed with the main result that 
g-Gaussians are indeed a good approximation of the numerical data if the parameter distance to the critical point and 
the number of iterations entering the sum satisfy a scaling condition that was derived in Q . In [l| the author also 
claims that a) Levy distributions could give an equally good fit to the data as g-Gaussians b) Levy statistics might 
have a better theoretical basis for this problem than g-statistics c) new types of distributions that he obtains by not 
neglecting transients could be universal. 

In this note we point out that the paper [H is misleading since most of the numerics performed in that paper 
operates in a parameter region that we explicitly excluded by the scaling condition derived in [3|. In the parameter 
region chosen by Grassberger, his statistics is insufhcient in the sense that much larger numbers of iterates would be 
needed to observe g-Gaussian distributions. Moreover we show that claims a) and b) are incorrect and that there is 
no theoretical or numerical basis for claim c). 

Let us use the same notation as in [ij . The object of study are sums Y of iterates Xi of the logistic map f{x) = 1 — 
with parameter a close to the critical point Qc = 1.4011551890920506... of period doubling accumulation. The sum 
consists of A'' iterates. One starts from an ensemble of uniform initial conditions and the first Nq iterates are ommitted: 

No+N 

Y= J2 (1) 

i=No + l 

The question investigated in [l|, 0, 0] is what probability distribution is to be expected for the random variable Y, 
since the ordinary Central Limit Theorem (CLT) is not valid close to the critical point due to strong correlations 
between the iterates. 

Before discussing the results of [H, we first point out a few formal errors in [l[. If a is slightly above the critical 
Qc then it is well-known that the attractor of the logistic map consists of n = 2*^ chaotic bands. In [l| it is stated (1 
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line below caption of Fig. 3) that k k [a — Oc)^/*, where 5 is the Feigenbaum constant. This statement is obviously 
wrong, the correct relation is 

There are a few further formal errors in 1]. In the caption of Fig. 4 in it is said that the figure shows numerical 
data for the probability density for various values of N and n. However, the actual data displayed in Fig. 4 seem to 
correspond to tuples of the form {n,N), i.e., the order of n and TV has been swapped. Another error is the fact that 
the absolute value is missing when the author refers to the z-logistic map fa,z{x) ~ a — \x\^ three lines after eq.(2). 

Let us now come to the actual content of the paper. In [3| we pointed out that the problem is more complex 
than the ordinary CLT, since two limits have to be performed simultaneously: a —^ ac and N — > oo. Simultaneous 
limits are standard knowledge in mathematics. In [3] we provided arguments that in order to obtain g-Gaussian limit 
distributions the simultaneous limit a ac and N ~^ oo must be performed in such a way that the scaling relation 

- 4*= (3) 

holds. Here k is again given by eq. and S = 4.6692011... is the Feigenbaum constant. In the notation used by 
Grassberger in [ll| our scaling condition is equivalent to 

N ~ (4) 

where again n denotes the number of chaotic bands. The abstract of [l| claims numerical inconsistency but the 
paper is misleading since most of the simulations in jlj ignore the above scaling condition ([3]) or but operate in 
a different parameter region (called 'peaked region' in Q). For example, for his results presented in Fig. 2 of [l| 
the author has chosen the distance \a — ad from the critical point Qc to be of the order 10~^^, basically fixed by 
his numerical precision. For this value our scaling relation gives k = 26.9 and hence ^ 4*^ ~ 2^'^ ^ 1.5 • 10^^ is 
required to see a g-Gaussian. On the contrary, the author performed his simulation in his Fig. 2 with the iteration 
numbers A^ = 256, 2048, 16384 and 131072, which are clearly insufficient to exhibit a g-Gaussian. Being that close 
to the critical point, much higher values of iteration numbers N are needed to obtain sufficient statistics to properly 
confirm or disconfirm our results presented in This is precisely the reason why in we chose larger distances 
from Oc for which the relevant A^- values are still reachable in a numerical experiment. 

In Fig. 4 and 5 of [l| the author investigates band splitting points and mainly looks at cases in the parameter space 
given by n > > VA or rt < < ^/N, again ignoring our scaling condition n y/N. Singular behavior of the density is 
observed simply because the scaling condition ([4]) is violated. Our present Fig. 1 shows how Fig. 4 of [H would have 
looked like had the scaling condition been satisfied. In this case one gets smooth curves that are well-approximated by 
q-Gaussians. In this sense, for example the band splitting point 4096 2048 investigated in [l|, which is the closest 
to the critical point, would yield the same g-Gaussian behavior had the iteration number A^ = 2^^ been used. 

Let us provide a few further comments. In [l[ some additional numerical experiments were performed testing the 
effects of different lengths Ao of omitted transients. The author emphasizes that in his opinion the condition No >> N 
is relevant. We have checked this claim in the relevant parameter region fixed by A^ ^ n^. As an example we have 
chosen one of the cases given in Fig. 2 of Q (namely, a — 1.401175 and A — 16384) and tested the effect of discarding 
transients of length A^o = 2048, 4096, 8192, 16384, 65536. The result is shown in Fig. 2. Apparently all curves fall onto 
each other, no matter whether Aq < A^ or Aq > A, and are well approximated by a g-Gaussian. Hence the condition 
Nq >> N advocated in [l| seems to be irrelevant in the scaling region, as long as Aq is sufficiently large (in 0, Q Aq 
was typically chosen to be 4096). 

In [l| it is also stated that Levy distributions, possibly motivated by the Levy-Gnedenko limit theorem, could give 
equally good fits. To test this claim, we tried to fit our data by Levy distributions as well. The result is shown 
in Fig. 3. The numerical data are well approximated by g-Gaussians, whereas Levy distributions give worse fits. 
More specifically, in a log- log plot. Levy distributions with 1 < a < 2 have an inflexion point which is by no means 
supported by the logistic-map data. If the parameter a of the Levy distribution is slightly increased, then the fit 
quality in the middle region is slightly improved but the tails become too pronounced to provide an acceptable fit. 
Hence the claim of [l| that Levy distributions might yield a better fit is incorrect. Besides this, the Levy-Gnedenko 
limit theorem holds for independent (or nearly so) random variables with infinite variance, whereas the iterates of the 
logistic map near to the critical point have strong correlations. Hence there is no theoretical justification to use Levy 
distributions in this problem. The conjecture in [l| that there might be a suitable ordering of the iterates into subsets 
that are almost independent lacks any theoretical proof or numerical justification. 

Finally, in Fig. 1 of Q, new types of distributions of Y are shown for the case Nq = 0, i.e. no transients are omitted, 
and claims are made at the end of the paper that these distributions including all the transients could be universal. 
As an argument for universality of transients at a = Cc, in [l| the work Q is cited. In however, only transient 
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FIG. 1: Investigation of the density of Y for various band splitting points with number of iterations A'' satisfying the scaling 
condition N ^ r? . 
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FIG. 2: Densities of Y obtained for various lengths A'o of omitted transients in the scaling region N -n? 



behaviour of iterates of the exact Feigenbaum fixed point function g is investigated, i.e. the map under consideration 
in [3| is the exact solution g of the Feigenbaum-Cvitanovic equation ag{g{x/a)) = g{x). However, universal behavior 
in our case would mean that different maps / with quadratic maximum would generate the same distribution of Y . 
Since different quadratic maps can have very different transient behavior, it seems highly unlikely that the sum Y of all 
these different transients would converge to a random variable that has a universal distribution, as claimed in For 
this, one would have to carefully estimate the speed of convergence of the iterates of / to the Feigenbaum fixed point 
function g under successive iteration and rescalation, which was not done in Hence there is no theoretical basis 
for the claim of that the observed transient distributions of sums are universal. Neither any numerical evidence of 
universality is provided in [l[ . 

Summarizing, the numerical experiments of [l[ were mainly performed in a different parameter region that was 
explicitly excluded by our scaling relation derived in Q . In the parameter region chosen by Grassberger the number 
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FIG. 3: Comparison of best fits obtained by using g-Gaussians (g — f .63) and Levy distributions (a — f .60), both in log-linear 
(a) and log- log plots (b). 



N of iterations is insufficient, fn the region fixed by our scaUng argument N ^ ri^ , q-Gaussians indeed provide good 
fits of the data, far better than the Levy distributions suggested in [1], which moreover do not have any theoretical 
justification for this problem involving strongly correlated random variables. Transient distributions investigated in 
m are unlikely to be universal. 

This work has been supported by TUBfTAK (Turkish Agency) under the Research Project number I04T148. C.T. 
acknowledges partial financial support from CNPq and Faperj (Brazilian Agencies). 



[1] P. Grassberger, Phys. Rev. E 79, 057201 (2009). 

[2] U. Tirnakli, C. Beck and C. Tsallis, Phys. Rev. E 75, 040106 (2007). 
[3] U. Tirnakli, C. Tsallis and C. Beck, Phys. Rev. E 79, 056209 (2009). 
[4] P. Grassberger and M. Scheunert, J. Stat. Phys. 26, 697 (1981) 



